Tile- Packing Tomography is NP-hard 



Marek Chrobak Christoph Diirr* Flavio Gumez^ Antoni Lozano"'- 

Nguyen Kim Thang§ 

December 22, 2010 



Abstract 

^v^j Discrete tomography deals with reconstructing finite spatial objects from their projections. The 

objects we study in this paper are called tilings or tile-packings, and they consist of a number of disjoint 
copies of a fixed tile, where a tile is defined as a connected set of grid points. A row projection specifies 
how many grid points are covered by tiles in a given row; column projections are defined analogously. 
For a fixed tile, is it possible to reconstruct its tilings from their projections in polynomial time? It is 

t-H known that the answer to this question is affirmative if the tile is a bar (its width or height is 1), while 

for some other types of tiles NP-hardness results have been shown in the literature. In this paper we 
present a complete solution to this question by showing that the problem remains NP-hard for all tiles 
other than bars. 

u 

oa 1 Introduction 

O 

Discrete tomography deals with reconstructing finite spatial objects from their low-dimensional projections, 
^vq Inverse problems of this nature arise naturally in medical computerized tomography, electron tomography, 

^ non-destructive quality control, timetable design and a number of other areas. This wide range of applications 

inspired significant theoretical interest in this topic and led to studies of computational complexity of various 
discrete tomography problems. For an extensive and detailed coverage of practical and theoretical aspects 
of this area, we refer readers to the book by Kuba and Herman, see jH [7] . 

In this paper we consider the problem of reconstructing a tile packing from its row and column projections. 
Formally, consider the integer grid of dimension m x n, consisting of all cells € [0, m) x [0, n). Every 
q>^ cell is adjacent to its neighbor cells 

(i-U), 

whichever of those are present in the grid. Alternatively, one can think of each (i,j) as a cell in an m x n 
matrix. In the paper, we will often use the matrix notation and terminology, using terms "row" and "column" , 
with rows numbered top-down and columns numbered from left to right, so that the upper-left cell is (0, 0). 

We define a tile to be any finite connected set T of grid cells. By "connected" we mean that for any 
two cells of T there is a path inside T between these cells, where any two consecutive cells on this path are 
adjacent. The width and height of T are defined in the obvious manner, as the dimensions of the smallest 
h x w rectangle containing T. If w = 1 or h = 1, then T is called a bar. By 

T+(i,j) = {(x + i,y+j):(x,y)eT} 
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we denote the translation of T by the vector T + (i,j) is called a (translated) copy of T, and 

is the position of this copy. Later in the paper, we will sometimes refer to T + as a "tile" - which is 

somewhat ambiguous but more intuitive than a "copy" 



s 




Figure 1: A tile packing of the 9 x 10 grid and its projections. By convention the leftmost upper cell is (0, 0). 

A tile packing of the m x n grid using T — or a T -packing, in short, if m and n are understood from 
context — is a disjoint partial covering of the grid with translated copies of T. Formally, a T-packing is 
defined by a set D of translation vectors such that all translated copies T + (i,j), for all G D, are 

contained in the m x n grid and are pairwise disjoint. We do not require the tiles to completely cover the 
grid. We stress here that what we call a T-packing here has been sometimes called a partial tiling in the 
literature; see, for example, [5]. An example of a T-packing is shown in Figure [TJ Without loss of generality, 
throughout the paper, we will be assuming that the tile T used in packing is in a canonical position in the 
upper-left corner of the grid, that is min {x : (x, y) G T} = min {y : (x, y) G T} = 0. 

To simplify notation, instead of counting how many grid points are covered by tiles in a given row (or 
column), we count how many tiles start in a given row (column), which is equivalent up to some base-change. 
So the row and column projections of a packing D are defined as a pair r G N m and s G N™ of vectors such 
that 

U = \{j:(i,j)€D}\ and Sj = \{i : G D}\. 

For example, consider tile T = {(0, 0), (0, 1), (1, 0), (1, 1), (1, 2), (2, 1), (2, 2)}. Figure [l] shows an example of 
a T-packing. This packing is 

D = {(0, 3), (2, 0), (2, 7), (3, 5), (4, 3), (6, 0), (6, 7)}. 

We study the problem of reconstructing tile packings from its horizontal and vertical projections. More 
formally, for any fixed tile T, the problem is defined as follows: 

Tile Packing Tomography Problem (TPTP(T)): The instance of TPTP(T) consists of vectors r G N m ,s G 
N". The objective is to decide if there is a T-packing D whose projections are r and s. 

This problem has been introduced in [2] and shown to be NP-hard for some particular tiles. In [T], the 
proof technique has been adapted to show NP-hardness for any rectangular tile, i.e. a tile that consists of 
all cells G [0, h) x [0, w) for some dimensions w, h > 2, 

On the positive side, the classical work of Ryser [TT] on projections of 0-1 matrices provides a character- 
ization of vectors that correspond to projections of T-packings for the special case when T is a single cell, 
and provides a simple polynomial-time algorithm for that case. The ideas from [TT] were extended in [31 E] 
to the case when T is a bar. In [TJ, polynomial-time algorithms were given for restricted special cases. The 
complexity status was unknown for all other tiles, and the current paper completes the picture by proving 
the following theorem. 

Theorem 1 Problem TPTP(T) is NV-complete for any tile T that is not a bar. 
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The general structure of our proof resembles those introduced in [5] and [3] , although the reductions we 
present are substantially more difficult, since the generality of our result means that we cannot take advantage 
of a specific shape of the tile, and that we need to base the construction on some generic properties shared 
by infinitely many tiles. Our techniques take advantage of Ryser's structure results for 0-1 matrices, in 
combination with some arguments based on linear algebra. 

After reviewing some background information in Section [2j we introduce the main idea of the reduction 
in Section |3j by formulating the overall framework of the reduction and conditions on T required for this 
reduction to be correct. Then, in Section [4] we show that all non-bar tiles satisfy these conditions. 

2 Main Tools 

In this section we briefly review two concepts that will play a crucial role in our proofs: affine independence 
and Ryser's theorem. 

Affine independence. Vectors Ui, V2, Vk € R n are called affinely independent if the unique solution of 
equations 

fe fc 

OLi = and = 

i=l j=l 

is ai = cti = ■■■ = cik = 0. It is easy to show that the following three conditions are equivalent: 

(ail) vi,V2-, Vk are affinely independent, 

(ai2) vi — vi, t?3 — Vi, Vk — v\ are linearly independent, 

(ai3) (ui, 1), (i>2, 1), {vk, 1) are linearly independent. 

We will refer to vectors Vi — «i, i = 2, 3, ...,k, in (ai2), as difference vectors. Condition (ai2) is useful in 
verifying affine independence. For example, (1, 1), (3, 4), (5, 5) are affinely independent because the difference 
vectors (3, 4) — (1, 1) = (2, 3) and (5, 5) — (1, 1) = (4, 4) are linearly independent. 

Condition (ai3) implies that if Vi,V2, •••,i>fc are affinely independent then for any vector v and constant 
(3, the equations 

fe fe 

otiVi = v and a,; = /3 

i=l i=l 

have a unique solution a.\, 012, etifc. 

For any vector v G K n and any set of indices ii>*2> € [0,n), define the (ii,«2: ■■.,%) -restriction 
of v to be the vector v' € M h that consists only of the coordinates i t , t = 1, b, of v. For example, the 
(0, 3, 4)-restriction of v = (4,3,1,0,7,9,5) is v' — (4,0,7). For any set of vectors vi,V2,— ) Vk G R n , to 
show that they are affinely independent it is sufficient to show that their (ii, 12, if,)-restrictions arc affinely 
independent, for some set of indices i\, %2, it- 

Ryser's theorem. Ryser [10] studied the structure of 0-1 matrices with given projections. We adapt his 
characterization of these matrices and express it in terms of tile packings. 

Fix a tile T and let / C [0, m) be a set of rows and J C [0, n) a set of columns. We say that a tile copy 
T + (i,j) belongs toIxJiii£l,j£j. Note that here we do not require inclusion of T + (i,j) in / x J. 

Define = max^i \D <~) (I x J)|, where the maximum is taken over all T-packings D of the m x n grid. 
Thus £/. j is the maximum number of copies of T that can belong to / x J in a T-packing without overlapping 
(and without any restriction on their projections). 

For a set / of rows, denote r(I) = J2iei r i- Analogously, let s( J) — J2jej s j> f° r a set J 01 columns. By 
/ = [0, m) — I and J = [0, n) — J we denote the complements of of sets / and set J, respectively. 
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Consider a T-packing D with projections r, s. Then we have 

r(I)-s(J) = \Dn (I x J) | - \Dn (J x J)\. 

By definition, \D PI (I x J)| < fy. Therefore we obtain the following lemma (inspired by |10j). 

Lemma 1 Let I be a set of rows and J a set of columns. If r{I) — s(J) = then every T-packing D with 
projections r, s satisfies \D PI (I X J)| = and \D D(I X J)\ = 0. 

3 General Proof Structure 

For each non-bar tile T, we show a polynomial-time reduction from the 3- Color Tomography Problem in- 
troduced in [6] and shown to be NP-hard in [4]. In that problem, an object to be reconstructed is a set 
of "atoms" (in our terminology, single cells) colored red (R), green (G) or blue (B). The instance contains 
separate projections for each color. The formal definition is this: 

3-Color Tomography Problem (3CTP): The instance consists of six vectors r R , r G , r B g N ,n , s R ,s G ,s B £ 
N". The objective is to decide whether there is an m x n matrix M with values from {R, G, B} such 
that, for each color c £ {R, G, B}, r% = \{y : M xy = c}\ for each x and s y = \{x : M xy = c}\ for each y. 

From now on, assume that T is some non-bar fixed tile of width w and height h. Let T be an instance 
of 3CTP for some m x n matrix specified by six projections r R , r G , r B , s R , s G , s B . We will map I into an 
instance J of TPTP(T) for an m! x n' grid with projections r, s, such that I has a matrix M with projections 
r R , r G , r B , s fl , s G , s B if and only if J has a T-packing with projections r, s. 

Without loss of generality we assume that for every color c we have ^ x r x = ^2 y s y , for every row x we 
have X) c r S = m ' aim -^ or ever y column y we have ^2 c s y = n. Otherwise, T is of course unfeasible, so we 
could take J to be any fixed unfeasible instance of TPTP(T). 

We now describe J . We will choose a grid of size m! x n' for m! — mk and n' = n£, where k and £ 
are positive integer constants to be specified later. We will use the term block for a k x I grid. We can 
partition our m! x n' grid into mn rectangles of dimension k x £, and we can think of each such rectangle as 
a translated block. The rectangle [xk, (x + l)k) x [y£, (y + 1)£) will be referred to as block (x,y). 

Next, we need to specify the projections r and s. We will describe these projections in a somewhat 
unusual way, by fixing three packings of a block denoted D R , D G , and D B (obviously, corresponding to 
the three colors), and then expressing r and s as linear combinations of these packings. More specifically, 
denoting by f c and s c the horizontal and vertical projections of packing D c , for each c € {i?, G, -B}, we define 

r xk +i = r * • f 1 and s v e +o =J2 s y'Sj C 1 ) 

c c 

for every i g [0, fc), j G [0,£), x € [0, m), and y G [0,n). The idea is that replacing each cell in a solution to 
the 3CTP instance r c , s c , for c = R,G, B, by its color-corresponding block gives a solution to the TPTP(T) 
instance r, s. 

To complete the description of the reduction, it still remains to define the three packings D R , D G , and 
D B . This will be done in the next section. In the remainder of this section we establish conditions that will 
guarantee correctness of our reduction. 

Our three packings will be designed to satisfy the following two requirements: 

Requirement 1: Vectors f R ,f G , r B are affinely independent and vectors s R ,s G , s B are affmely independent. 
Note that, by property (ai3), this implies that for any vector v there is at most one possible way to 
represent it in a form v = nnf R + n<jf G + nsf 3 , where ur + ng + tib = n. Naturally, an analogous 
statement holds for column projections. 

Requirement 2: In any packing D of J with projections r, s, the restriction of D to each block of the grid 
has projections equal to r c , s c , for some c £ {R, G, B}. 
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Lemma 2 Assume that the three packings D R , D G , D B satisfy Requirements 1 and 2. Then! has a solution 
if and only if J has a solution. 

Proof: Let M £ {R, G, B} mxn be a solution to 1. We transform M into the following packing D 
for the m! x n' grid: 

D = U U (D M " + (xk,y£)). 

i6[0,m) y€[0,n) 

In other words, if M xy = c then block (x, y) of the m' x n' grid contains a copy of D c . By simple inspection, 
the projections of D are indeed equal to the vectors r and s in (JTJ. 

(•<=) For the converse, suppose that there is a packing _D with projections r, s. By Requirement 2, every 
block of the m! x n' grid has projections f c and s c , for some c G {ii, G, B}. We then associate this block 
with color c. We can thus define a matrix M £ {R,G, B} mxn such that M xy = c if block (x,y) of 13 is 
associated with color c. 

We now need to show that M is a solution for X. To this end, fix some arbitrary < x < m and consider 
vector 

v = (r xk ,r xk+1 , ...,r( x+ i)fc_i), 

which is the projection of the "row" of all blocks (x, y), for all y. By the construction ([!]), v can be written 
as 

v = nnf R + nof G + jibt b , 

where ur — r x , uq — r x , and ns = r x . Now, using Requirement 1, we obtain that this representation is 
unique under the assumption that nj^ + uq + ng — n. We can thus conclude that the projection of row x 
of M is correct, that is \{y : M xy = c}| — r x for all c. By the same argument, column projections of M are 
correct as well, completing the proof. □ 

In summary, to complete the proof for the given tile T, we need to do this: (i) define a rectangular 
k x t block with three packings D R , D G , D B , (ii) show that the row projections of D R , D G , D B and the 
corresponding column projections are affinely independent (Requirement 1), and (iii) show that in any 
solution to J ', each block (x, y) has projections equal to those of one of D c , for some c (Requirement 2). We 
show the construction of such block packings in the next section. 



4 Construction of Block Packings 

As in the previous section, T is a fixed (but arbitrary) non-bar tile. We call (i,j) a conflicting vector if T 
and T + (i,j) overlap, that is Tfl (T + ^ 0- Obviously, the vectors (i,j) and are either both 

conflicting or both non-conflicting. Since T is not a bar, it has a conflicting vector with i, j ^ 0. To see 
this, observe that since T is connected and not a bar, it contains two cells (i,j) and {il with \i — i'\ = 1 
and \j — j'\ = 1, so one of (1, 1) or (—1, 1) is a conflicting vector. 

For the construction of the proof, fix a conflicting vector of T that maximizes the Li norm under the 
constraint that none of the coordinates is 0. We denote this vector by (— p, q) and assume without loss of 
generality p, q > 0, for otherwise we can flip T horizontally or vertically and give the proof for the resulting 
tile. 

So any vector with i, j ^ and \i\ + \ j\ > \p\ + \q\ is not conflicting. Let a be the smallest positive 
integer such that (ap, 0) is not a conflicting vector. Similarly let b be the smallest positive integer such that 
(0, bq) is not a conflicting vector. Without loss of generality we assume that a ^ b, since otherwise we can 
exchange the roles of columns and rows in the proof. 

We now divide the proof into four cases, and for each of them we show that Requirements 1 and 2 are 
satisfied. 
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Figure 2: For each case of the proof, the figure shows three packings (left-hand side) and an example of a tile 
satisfying this case (right-hand side). Dotted vectors are non-conflicting, either by maximality of (— p, q) or 
by the case condition. In the third case, all vectors that are not shown (to avoid clutter) are non-conflicting. 



4.1 Case a = 1 and b = 1 

In this case, we use the following three packings: 

D R = {(p,0),(p,q)}, 
D G = {(0,q),(p,q)}, 
D B = {(p,q)}. 

The values of k and I are chosen to be the smallest integers for which these three packings are contained in 
the k x I grid. 

The packings are depicted on Figure [2j The squares represent possible positions for tiles. Two positions 
are connected with a solid line if the difference of the positions is a conflicting vector. That means that no 
packing can contain simultaneously a tile in both positions. Dotted lines indicate non-conflicting vectors, i.e. 
they connect pairs of compatible positions. We show lines only between position pairs relevant to the proof. 
In the figure we mark with letter c the tile positions of D c , for c 6 {R, G, B}. For illustration, for each case, 
on the right-hand side of the figure we show a tile satisfying the case conditions. Again, solid vectors are 
conflicting and dotted vectors are non-conflicting. 

We first verify Requirement 1. The (O,p)-restrictions of f R ,f ,f B are, respectively, (0,2), (1,1) and 
(0, 1), and the (0, (^-restrictions of s R , s , s B are, respectively, (1, 1), (0, 2) and (0, 1). For both the row and 
column projections, routine calculations show that their restrictions are affinely independent. 

We now focus on Requirement 2. Let r, s be the projections obtained by the reduction, and consider 
a packing D with these projections. We use Lemma [T] with / being the set of all row indices that are p 
modulo k, and J being the set of all column indices that are q modulo £. By inspecting the definition of the 
projections we have 

r(I) = r R + mn = s R + ran = s(J) + mn, 

x£[0,m) ye[o.n) 
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so r(I) — s(J) — mn, which is |7 X J|. In consequence Lemma [T] applies. Therefore every block in D contains 
a tile at position (p, q) and none at position (0, 0). The remaining possible positions for tiles are (p, 0), (0, q), 
but both cannot be occupied in the same block. This shows that every block of D is one of the packings 
D R ,D G ,D B . 

4.2 Case a = 1 and b > 2 

In this case, we use the following three packings: 

D R = {(2p,0),(2p,bq),(p,0)} 
D G = {{2p,Q), (2p,bq),(0,q)} 
D B = {(2p,0),(2p,bq)} 

The values of k and £ are chosen to be the smallest integers for which these three packings are contained in 
the k x I grid. 

The (0,p, 2p)-restrictions of f R ,f G ,f B are linearly independent vectors (0,1,2), (1,0,2), (0,0,2); there- 
fore f R ,f G ,f B are affmely independent. By a similar argument, we obtain that the corresponding column 
projection vectors s R ,s G ,s B are affmely independent as well. Thus, Requirement 1 holds. 

Now we verify Requirement 2. Row 2p of a block can contain at most 2 tiles. Let I be the set of rows i 
with i mod k = 2p. By inspecting the definition of the projections we have r(J) = 2mn, so every block in a 
solution D must contain exactly 2 tiles in row 2p, and they are at positions (2p, 0), (2p, bq). 

Now let J be the set of all columns j with j mod £ = bq. We only have s( J) = mn, so in every block of 
D the positions (0, bq), (p,bq) are empty. The tile at (2p, 0) forces position (p,q) to be empty. By the case 
assumption that a = 1, there is no conflict between (2p, 0) and {p, 0). By maximality of (~p,q) there is no 
conflict between positions (2p, 0) and (0, q), or between (2p,bq) and (p, 0), or (2p,bq) and (0, g). 

That leaves only 3 positions where the block packings can differ, namely (0, 0), (p, 0) and (0, q). Let d be 
the number of blocks in D with a tile in (0, 0). Similarly let e be the number of blocks in D with a tile in 
(0, q). Now we use the fact that in the original instance we had ^ r G = ^ s G ; let Ag denote this quantity. 
This time, let /' be the set of rows i with i mod k = 0, and J' the set of columns j with j mod £ — q. Then, 
by the definition of d and e, we have r(J') = d + e, and by the definition of the chosen three packings, we 
have r(I') — Ac- Similarly, we have r(J') = e and r(J') — Ag, which shows d = 0. Therefore every block 
packing in D is one of D R , D G or D B . 

4.3 Case a = 2, 6 = 2 and Vector (p,q) Conflicting 

In this case we assume a = 2, b = 2 and that the vector (p, g) is conflicting. Since (— p, q) is conflicting as 
well, this makes the construction very symmetric. The three packings used in this case are: 

D R = {(0,2q),(p,0),(2p, 2q)}, 
D G = {(Q,q),(2 Pl Q),(2p, 2q)}, 
D B ={(0,0),(0,2 9 ),(2p,0),(2p,2?)}. 

Again, the values of k and £ are chosen to be the smallest integers for which these three packings are contained 
in the k x £ grid. The construction is illustrated in Figure [2j The idea behind this construction is similar to 
the reduction used in [2] to show NP-hardness of the packing problem for the 2x2 square tile. 

The (0,p, 2p)-restrictions of f R ,f G \f B are (1,1,1), (1,0,2), (2,0,2); therefore f R ,f G 1 f B are affmely 
independent. By symmetry, the same holds for s R ,s G ,s B . 

Now we verify Requirement 2. Let D be a packing with projections r, s. Due to conflicts, there can 
be at most 2 tiles in rows p, 2p of a block. Let I = {i : i mod k £ {p, 2p}}. Since r(J) = 2mn and there 
are mn blocks in D, every block of D must contain exactly 2 tiles in rows p, 2p. By symmetry, the same 
holds for columns q, 2q. There are only 4 packings that satisfy these constraints and avoid conflicts, namely 
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D R , D G , D B and packing D A = {(0, 2q), (2p, 0), (2p, 2q)}. Let Xr, Xg, Xb, Xa be the respective numbers of 
these different block packings in D. Since D has row projection r, by expressing the total number of tiles in 
two different ways, we have 

3 X! r ^ + 3 X^ +4 X^ =3X r + 3X g +AXb + 3X a . (2) 
From the assumption that ^2 C ^2 X r x — mn, we have also 

J2 r * + H r x + H r * = ™ n = X R + X G + X B + X A . (3) 

Now let I' = {i : i mod k = p}. Then ^ x r R — r(I') = Xr. Similarly, for J' = {j : jmod I = q} we obtain 
J2x r x — J2 y s y = S (J') = ^g- After subtracting these two equations from |2| and (|3j), we are left with two 
equations 

J2^ = X B +X A , 4j]rf =4A b + 3A j4 , 

from which we conclude Xa = 0. This verifies Requirement 2. 

4.4 The Remaining Case 

Assume now that none of the previous cases holds. Since a < 6, this means that either 

(a) a > 2 and b > 3, or 

(b) a — 2, b = 2 and vector (p, q) is not conflicting. 

We claim that the vector (p, (b — l)q) is not conflicting. Indeed, in case (b) above, it follows by case 

assumption and in case (a), it follows from the maximality of (—p, q). Therefore, in any block of a T-packing 
both positions (0, q) and (p, bq) could contain a tile. 
We use the following three packings (see Figure [2]) : 

D R = {(p,0),(p,bq)}, 

D G = {(0,q),{p,bq)}, 
D B = {(p,bq)}. 

Again, the values of k and £ are chosen to be the smallest integers for which these three packings are contained 
in the k x £ grid. 

The (0, ^-restrictions of f R ,f ,f B are (0,2), (1,1), (0,1), and their difference vectors (1,-1), (0,-1) are 
linearly independent. Therefore r R ,f G ,f B are affinely independent. The (0, q, 6(?)-restrictions of s R ,s G ,s B 
are linearly independent vectors (1,0,1), (0,1,1), (0,0,1); therefore s R ,s G ,s B are affinely independent. 
Thus, Requirement 1 holds. 

Now we verify Requirement 2. The proof is similar to those of previous cases, except that now we have 
more candidate packings to consider. Fix some T-packing D with projections r, s. First, using the same 
arguments as in the previous case, we observe that every block of D must contain exactly one tile in column 
bq, that is either in location (0,bq) or (p,bq). Taking conflicts into accounts, straightforward case analysis 
produces nine possible packings, including two pairs of "equivalent" packings with identical projections. We 
now introduce notation for the numbers of these packings: 

A /v. = the number of packings D R — {(p, 0), (p, bq)}, 

Xg = the number of packings D G = {(0, q), (p, bq)} or {(p, q), (0, bq)}, 

Xb = the number of packings D B = {(p, bq}), 
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Xa = the number of packings {(0, 0), (p, bq)} or {(p, 0), (0, bq)}, 

Xc = the number of packings {(0, 0), (p, q), (0, bq)}, 

Xd = the number of packings {(0,bq)}, 

Xe = the number of packings {(0, 0), (0, bq)}. 

Let I be the set of all rows i with i mod k — and J be the set of all columns j with j mod I = q. 
Then, by inspecting the projection definitions, we have r(I) — ^2 x r x = ^2 y s^ = s(J). Since r and s arc 
the projections of D, we also have r(I) = Xq + Xa + 2Ac + Xe + 2Xe and s(J) = Xc + Xc- This shows 
Xa = Xc = Xe = Xe = 0, completing the analysis of this case and the proof of NP- hardness. 
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